Flux surface shaping effects on tokamak edge turbulence and flows 
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Shaping of magnetic flux surfaces is found to have a strong impact on turbulence and transport in 
tokamak edge plasmas. A series of axisymmetric equilibria with varying elongation and triangularity, 
and a divertor configuration are implemented into a computational gyrofluid turbulence model. The 
mechanisms of shaping effects on turbulence and flows are identified. Transport is mainly reduced 
by local magnetic shearing and an enhancement of zonal shear flows induced by elongation and 
X-point shaping. 
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I. INTRODUCTION 



It can be regarded as one of the most intriguing re- 
sults of half a century in physics of hot magnetized plas- 
mas, that in macroscopically stable equilibria the con- 
finement is nevertheless determined by micro-scale insta- 
bilities and fluctuations, with transport of energy and 
particles across magnetic flux surfaces being dominated 
by ubiquitous turbulence on small (gyro radius) scales 
and low (drift) frequencies [1 S II llB IE ■ This tur- 
bulence is regulated by the formation of mesoscopic zonal 
structures out of the turbulent flows |sl 1^. ITol Hit Il2t Il3| . 

Criteria for design of new magnetized plasma experi- 
ments for fusion research like ITER are primarily based 
on empirical scaling laws for the energy confinement time 
[i~H |. First-principle based transport models still require 
reference to experimental scalings, and are validated 
mainly only for core plasmas excepting the pedestal re- 
gion |5j. An improved understanding of the underlying 
turbulent dynamics and of its relation to design param- 
eters, like the shape of the confining magnetic field, will 
facilitate the development of advanced magnetic confine- 
ment experiments. 

The plasma shape of a tokamak enters into confine- 
ment and transport modelling through parameters spec- 
ifying a vertical elongation n > 1 and an outboard side 
triangularity 5 > that describe the deviation from a 
simple circular torus. The design criteria of the ITER 
plasma shape with k = 1.7 and S = 0.33 (at 95% flux 
surface) are based on a conservative regime that is well 
established in present experiments More extreme 

variations of flux surface shaping were investigated with 
the experiment TCV (Tokamak a Configuration Vari- 
able) which is able to achieve configurations up to n = 3 
and 6 = ±0.5 [l5|. Experimental evidence suggests that 
large elongation is always beneficial for local and global 
energy confinement, whereas the role of triangularity de- 
pends also on the value of the plasma pressure gradient. 



FIG. f : Poloidal cross section of a tokamak plasma (ASDEX 
Upgrade) showing the elongation and triangularity of closed 
flux surfaces and a lower X-point on the separatrix. A sec- 
tion of the flux tube for drift wave turbulence computations 
winding around the torus along field lines is here locally in- 
tersecting the outer midplane of the plasma edge. 

An important task for computational plasma physics 
is to determine the effects of flux surface shaping on tur- 
bulence and transport by first-principle direct numerical 
simulation. Shaping effects are complicated even for ap- 
parently simple circular equilibria [lg , and efforts to per- 
form modelling for edge plasmas most often incorporate 
technically complicated prescriptions in an inflexible way 
|17| . An advantageous way to explore shaping is through 
the use of a covariant metric in a model with Cartesian 

topology mm, 

but even then it was found that the use 
of experimental geometry, while most realistic for mod- 
elling purposes, may not be helpful to theoretical under- 
standing, and hence standard "kappa-delta" models were 
used |2J|. In attempt to understand the effects of stel- 
larator geometry, a set of analytical forms were used for 
the metric, varying each in turn to provide a small set of 
steps from tokamak to stellarator geometry [2l|. Forms 
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like this were also used to characterise the X- point effects, 
again making simple control tests tractable |22j] . 

In the following, we report on the results of a more sys- 
tematic series of "kappa-delta" tokamak configurations 
with varying elongation and triangularity, and an ac- 
tual ASDEX Upgrade [23| divertor configuration with 
the equilibrium code HELENA [2^|. These are imple- 
mented into the gyrofluid code GEM ([25J, which is up- 
dated from Ref. [2(| ) for computations of edge turbulence 
and flows. Specific effects of flux surface shaping on drift- 
Alfven (DALF) and ion temperature gradient (ITG) tur- 
bulence in the tokamak edge plasma are analysed and 
compared, and the mechanisms of the coupling between 
geometry and turbulence are identified. The emphasis is 
on physical mechanism, which is complimentary to and 
should also be helpful to the understanding of the more 
experimentally oriented modelling efforts. 

It is found that turbulent transport fluxes are reduced 
by increasing flux surface elongation n. The influence 
of triangularity S is much weaker and generally also de- 
pends on elongation and on a pressure shift. The geodesic 
acoustic oscillation modes (GAM) and the overall turbu- 
lence frequency spectra depend on flux surface geometry. 
An enhancement of zonal flow amplitude in the spectra 
by elongation and X-point shaping is found to be weak 
in the DALF case and more pronounced in the ITG case. 
The reduction of transport observed in our simulations 
in shaped tokamak geometry is mainly a result of local 
and global magnetic shear. 

We begin in Section II with a brief review and discus- 
sion of the role of sheared flows in tokamak turbulence in 
context of the toroidal flux surface geometry. In section 
III the gyrofluid electromagnetic model GEM for edge 
turbulence computations is presented, and the flux tube 
representation of tokamak geometry within a GEM tur- 
bulence code is discussed in Section IV. Transport scal- 
ings by elongation and triangularity, that are obtained 
with this model, are presented in Section V. Resulting 
frequency and wave number spectra in tokamak divertor 
geometry are discussed in Section VI. Conclusions and 
implications of the results are presented in Section VII. 



II. THE ROLE OF SHEARED FLOWS IN 
TOKAMAK TURBULENCE 

One of the main requirements for successful perfor- 
mance of fusion experiments is the prospect of operation 
in a high confinement H-mode The forma- 

tion of a characteristical transport barrier in magnetically 
confined toroidal H-mode plasmas is closely related to 
the presence of a radially varying electric field near the 
last closed flux surface 27, 29]. A radial electric field, 
E r , gives rise to a perpendicular E-cross-B (ExB) flow 
through the relation v^; — (c/B 2 )ExB. If E r varies with 
r, then this flow is sheared, and this flow shear is what 
is invoked to explain the transport barrier [271 l29l l3fj| . 
The flow shear (vorticity) is expected to suppress tur- 



bulent transport either by shearing the eddies apart |3lJ 
or by simply tilting them so that their Reynolds-stress 
interaction with the vorticity layer tends towards energy 
transfer out of the eddies and into the flow Neo- 
classical equilibrium effects can also produce these shear 
layers H^HH]. As opposed to these mean ExB flows, the 
fluctuating zonal flows in direct interaction with the tur- 
bulence have also been observed in tokamak edge regions 

M 

The transition between states of low (L) plasma con- 
finement to the high (H) confinement regime featuring 
such an edge transport barrier is up to now not satis- 
factorily described by any of the existing first-principle 
theories |3q |. The notion that generation of zonal flows 
by direct drive of the turbulence can act to self-regulate 
by suppressing the driving turbulent vortices has lead to 
the development of predator-prey type bifurcation mod- 
els that are able to describe specific characteristics of the 
L-H transition (see e.g. Early computations based 

on the resistive-g [3£| and collisional drift wave [39] mod- 
els found self-generated sheared flows which could be ar- 
gued as an L-H transition trigger. However, more recent 
models in more comprehensive treatment of the toroidal 
geometry and also the full range of scales of motion do 
not find this |20| . The phenomenology of these self con- 
sistent zonal flows receives wide interest 

mm nil. but 

they serve as a moderator of the turbulence rather than 
a suppression mechanism, as the energetic interactions 
with not only the pressure but also magnetic parts of the 
equilibrium tend to place each part of the overall dynami- 
cal system, which naturally includes the flows themselves 

m- 

After the initial discovery of the H-mode in the divertor 
experiment ASDEX j2(| it has been found in numerous 
large tokamaks around the world that the presence of 
a divertor (originally designed to improve the impurity 
exhaust) significantly enhances the prospect of reaching 
this state [41|. However, up to now both the transition 
models as well as the numerical simulations have relied 
mainly on simplified representations of the plasma ge- 
ometry. In our present work we investigate the effects of 
realistic flux surface shaping on tokamak edge turbulence 
using a gyrofluid model recently corrected for energy con- 
servation properties 0, |2j|, thereby updating previous 
gyrofluid results [2(j. This provides a counterweight to 
existing efforts based on collisional fluid equations and 
associated numerical schemes which cannot treat the ion 
gyroradius scales and therefore overestimate linear phe- 
nomena based in the larger scales [TtI l42| . The ion gy- 
roradius necessarily enters because the dynamics is sen- 
sitive to phenomena in the range of the spectrum with 
perpendicular wavenumbers only slightly less than the 
inverse of the drift scale (p s , also called "ion sound gy- 
roradius"), and with comparable ion and electron tem- 
peratures p s is commensurate with the ion gyroradius pi 
[f|. A schematic view of a shaped plasma cross section 
including a flux-tube section from turbulence computa- 
tions is shown in Fig. ^ 
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Transport in a fusion edge plasma is dominated by 
turbulent low-frequency drift wave motion that causes a 
fluid-like convection through ExB vortices in a plane per- 
pendicular to the magnetic field direction, acting in con- 
cert with parallel coupling dynamics tending towards dis- 
sipation [43l 1441143 . In general this is both nonlinear and 
electromagnetic [(J Q . Toroidal compressibility of zonal 
(flux surface averaged) ExB flows opens channels in the 
nonlinear dynamics of the vortex-flow interaction that 
strongly affect the flow energetics 0, 0] . The origin of 
the compression is the geodesic curvature of the magnetic 
field lines. The flow is thereby coupled with poloidally 
asymmetric pressure sidebands which are consumed by 
the turbulence and global Alfven or parallel flow dynam- 
ics. This energetics was given in Ref. |47j and the various 
coupling effects throughout the flow/equilibrium system 
were analysed in Refs. pfll |4S[ . and the sensitivity of 
the turbulence to the geodesic energy transfer effect was 
explored using an additional scaling parameter for the 
strength of its coupling [43 ■ 

This geodesic transfer mechanism represents a restor- 
ing loss channel for the zonal flows, ultimately placing 
them in statistical equilibrium with the turbulence, with 
the Reynolds stress (spin-up) mechanism continuing to 
operate. Since the latter acts on all frequencies, the 
geodesic acoustic oscillation (GAM) itself need not 
be present for the geodesic transfer mechanism to oper- 
ate. The net result of this is that turbulent transport 
in 3D toroidal edge computations is found to be reduced 
but not completely suppressed by the self-generated zonal 
flows. This geodesic transfer effect was first studied 
for the case of a simple circular toroidal magnetic field 
[40L I47I l48j , In the following we give both numerical evi- 
dence and a physical picture of how plasma shaping and 
the presence of a divertor via its modification of plasma 
geometry can influence the outcome of the gyrofluid tur- 
bulence/flow interaction. 



III. THE GYROFLUID ELECTROMAGNETIC 
MODEL GEM 

The dynamical character of cross-field ExB drift wave 
turbulence in the edge region of a tokamak plasma is gov- 
erned by electromagnetic and dissipative effects in the 
parallel response. The most basic drift-Alfven (DALF) 
model to capture the drift wave dynamics includes non- 
linear evolution equations of three fluctuating fields: the 
electrostatic potential </>, electromagnetic potential A\ | , 
and density n. The tokamak edge usually features a 
more or less pronounced density pedestal, and the dom- 
inant contribution to the free energy drive to the turbu- 
lence by the inhomogeneous pressure background is thus 
due to the density gradient. On the other hand, a steep 
enough ion temperature gradient (ITG) in the edge does 
not only change the turbulent transport quantitatively, 
but adds new interchange physics into the dynamics. In 
addition, more field quantities have to be treated: par- 



allel and perpendicular temperatures Tn and T± and the 
associated parallel heat fluxes, for a total of six moment 
variables for each species. Finite Larmor radius effects 
introduced by warm ions require a gyrofluid description 
of the turbulence equations. Both the resistive DALF 
and the ITG models can be covered by using the six- 
moment electromagnetic gyrofluid model GEM [2{| , but 
we will refer for the DALF model to its more economical 
two- moment version 

The gyrofluid model is based upon a moment approxi- 
mation of the underlying gyrokinetic equation. The first 
complete six-moment formulation was given for slab ge- 
ometry |51|, and later extended to incorporate toroidal 
effects [52| using a ballooning-based form of flux surface 
geometry j5^|. Electromagnetic induction and electron 
collisionality were then included to form a more general 
gyrofluid for edge turbulence |20j ] , with the geometry cor- 
respondingly replaced by the version from the edge tur- 
bulence work, which does not make ballooning assump- 
tions and in particular represents slab and toroidal eigen- 
mode types equally well and does not require radial pe- 
riodicity |l9j| . Energy conservation considerations were 
solidified first for the two-moment version [(|, and most 
recently for the six-moment version [25J. This latter ref- 
erence currently defines the GEM model, energetics, and 
also the numerical techniques. The present paper uses 
the six-moment equation set for each species (Eqs. 99- 
104), induction and polarisation equations (Eqs. 89,92), 
the zero gyroradius limit for electrons, the local model 
with separate linear drive terms and Dirichlet bound- 
ary conditions (Sec. IV A), Pade forms for gyroaveraging 
operators (Sec. IV B), and the dissipation model and 
numerical scheme (Sees. IV C-E), all with equation and 
Section numbers from Ref. [2f|. This is the standard 
local form for the GEM code. 

The flux surface geometry of a tokamak enters into 
the gyrofluid equations via the curvilinear generalisation 
of differentiation operators and via inhomogeneity of the 
magnetic field strength B. The different scales of equilib- 
rium and fluctuations parallel and perpendicular to the 
magnetic field motivate the use of field aligned flux co- 
ordinates. The differential operators in the field aligned 
frame are the parallel gradient 

V|| = (1/B)(B + Bl)-V, (1) 

with magnetic field disturbances Bj^ = (— 1/_B)B x VA11 
as additional nonlinearities, the perpendicular Laplacian 

Vi = V- [{-1/B 2 )B x (B x V)], (2) 

and the curvature operator 

/C = V-[(c/B 2 )BxV)]. (3) 

The dynamical character of the system is further de- 
termined by a set of parameters characterising the rel- 
ative role of dissipative, inertial and electromagnetic ef- 
fects in addition to the driving by gradients of density 
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and temperature. In particular, we specify collision- 
ality C — 0.51e(i/ e c s /Lj_)(m e /M i ), magnetic induction 
(3 = e(AiTp e / B 2 ), electron inertia fi — e(m e /Mi) and ion 
inertia e = (qR/ L±) 2 , where L± is the background gra- 
dient scale length, and c s = yjT e /Mi is the sound speed. 
We use a standard set of edge parameters given as C = 5, 
/3=l,A = 5ande = 18350. 

The perpendicular scale length for the DALF model is 
set as L± = L n to the density gradient length, and for 
the ITG model as L± = Lt% — 0.5L„ = 0.5Lx e to the ion 
temperature gradient length, so that rji = L n /LTi = 2. 
The ITG model has otherwise identical parameters. Ac- 
curate experimental measurement of ion temperature 
profiles in the tokamak pedestal region are only recently 
being advanced and have previously been plagued by 
large uncertainties. The present experimental knowledge 
of edge parameters can thus be seen commensurate with 
both the DALF and ITG model cases used here. 

The computational domain is set to 64 x 256 nodes in 
units of the drift scale p s = (c/eB)\/T e Mi for (x,y) and 
16 nodes in one field line connection length (2?tqR) in 

— 7T < Z < 7T. 

The dimensions are chosen to appropriately account 
for statistical isotropy in small scales in both perpendic- 
ular directions with satisfactory spectral overlap, and for 
an extended box size in the electron drift direction y to 
achieve convergence. A box resolution down to p s guar- 
antees inclusion of all scales necessary for the nonlinear 
dynamics that are esssential for the drift wave turbulence 
characteristics. 



IV. FLUX TUBE REPRESENTATION OF 
TOKAMAK GEOMETRY 

Tokamak equilibria are computed by solving the Grad- 
Shafranov / Lust-Schluter equation with the code HE- 
LENA 24] , implemented as described in Ref. [2(| • A set 
of nested flux surfaces in straight field line Hamada coor- 
dinates (V, 8, C) is obtained by specification of given ra- 
dial profiles of pressure and rotational transform, and of 
the shape of the bounding last closed flux surface. These 
Hamada coordinates are then transformed into a field 
aligned system and re-scaled into local flux tube coor- 
dinates (x,y,z) |18|. Global consistency in the parallel 
boundary condition for one connection length is main- 
tained 54]. A transformation of the coordinate y, which 
signifies the electron diamagnetic drift direction perpen- 
dicular to the magnetic field within flux surfaces, is ap- 
plied in order to avoid grid deformation by local magnetic 
shear [l9|. Otherwise, grid cells are sheared strongly in y 
direction particularly near the X-point region, with ma- 
lign consequences on nonlinear dynamics, especially in 
the vorticity, that lead to a violation of the basic drift 
wave character and overestimate linear MHD-like dynam- 
ics. The differential operators are then expressed in terms 
of the flux tube coordinates: The curvature operator be- 
comes 

JC = JC e (z)V x +K?(z)V y , (4) 
the perpendicular Laplacian in flute mode ordering is 
Vi - g xx (z)d xx + 2g x y(z)d xy + g yy (z)d yy , (5) 
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FIG. 2: Metric element g xx , g xy (unshifted), normal curva- 
ture K? , geodesic curvature !C X , magnetic field strength B and 
local magnetic shear S in a tokamak with elongation k = 1 
and triangularity 5 — (dashed lines), compared to a config- 
uration with k — 2 and 5 — 0.4 (thin solid lines) and to an 
actual ASDEX Upgrade configuration (bold solid lines). 



and the parallel derivative is 

V|, = b*(z)d z , (6) 

noting that the factor of B 2 in p 2 also depends on z. 
Some metric coefficients g lJ that were obtained for elon- 
gation k = 1 and K = 2 with triangularity 5 = and 5 = 
0.4 are shown in Fig. [21 Increasing elongation k specifi- 
cally rises the local magnetic shear S = V|| (g xy /g xx ) and 
reduces JC X both in the upper and lower regions of the 
torus that correspond to flux tube coordinates z = ±7r/2. 
In the following, the configuration with k = 1 and 5 = 
will be refered to as a "simple circular torus" (SCT). 

Local and global magnetic shear are in general known 
to have a damping influence on tokamak edge turbu- 
lence [2^ , whereas geodesic curvature acting through K. x 
upon the axisymmetric component (k y = components, 
or "modes") maintains the coupling for a loss channel 
from zonal flow energy eventually to turbulent vortices 
pTf. Both mechanisms help to reduce turbulent trans- 
port. Normal curvature JC V on the other hand strength- 
ens primarily the interchange forcing of the turbulence 

{ky + 0). 
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V. TRANSPORT SCALING BY ELONGATION 
AND TRIANGULARITY 

A series of tokamak equilibria is constructed for elonga- 
tion 1 < k < 2 and triangularity < <5 < 0.4 with equal 
profiles of pressure and rotational transform. The flux 
tube is chosen at a radial position r = \/V/V = 0.95, 
where Vq is the volume enclosed by the last closed flux 
surface. The background density profile and rotational 
transform for the GEM turbulence computations are lin- 
earised within the bounds of the radial computational do- 
main. We restrict our simulations for now to the closed 
field line region lying a few tens of ion gyroradii inside 
from the separatrix, thus avoiding complications that oc- 
cur by a divergent metric and associated grid deforma- 
tion when the last closed flux surface is approached in the 
vicinity of an X-point. Ultimately, the goal of edge tur- 
bulence simulations will be to combine sufficiently well 
resolved nonlocal drift wave computations with a repre- 
sentation of the realistic field line geometry crossing the 
separatrix to the bounded scrape-off layer region. 

When the plasma shape is thus varied, we have a par- 
ticular interest in the effects on fluctuation time and spa- 
tial scales, on the radial variation of flux surface aver- 
aged (zonal) flows, and on turbulent transport. We first 
apply the ensemble of shaped tokamak equilibria to the 
DALF model. The cross-field turbulent electron trans- 
port is determined by the flux surface average of the 
radial E x B convection of the fluctuating density by 
F e = {v x n e ) yz , given in standard gyro-Bohm normalisa- 
tion to n e c s (p s /L±) 2 . The transport flux F e {t) is fluc- 
tuating in time, and for specifying a quantitative value 
to it the time average over a sufficiently long window, 
that covers all relevant frequency scales, is taken during 
the final phase of simulations after all initial linear tran- 
sients and spin-up of flows and oscillations have reached 
saturation. 

As shown in Fig. [31 we find that F e is reduced to around 
a third when elongation is increased from a circular cross 
section to k — 2. Scaling F e against triangularity 5, 
we find in Fig. ^ that for n = 1 the transport flux is 
independent from 6 within the error bars given by the 
deviation of the fluctuating data from its time average. 
For higher elongation k = 2 the transport F e is increased 
from the normalised value of 0.20 to 0.28 by increasing 
5 from to 0.4. A good fit to the transport data in the 
DALF model within fluctuation error bars is obtained by 

F e DALF (n,S)^F h0 -^ 2 - 15+1A65 (7) 

where Fi o = F(k=1, (5 = 0). This reduction of transport 
by elongation can be related to three different causes: 

(1) Increase of local and global magnetic shear (mag- 
netic shear damping), (2) reduction of curvature drive 
(interchange coupling), and (3) an enhancement of 
sheared zonal flows (geodesic transfer). 
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FIG. 3: DALF: Dependence of turbulent electron transport 
F e (in gyro-Bohm units) on elongation k for triangularity S = 
0, 0.2 and 0.4 (from bottom to top). 
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FIG. 4: DALF: Dependence of turbulent electron transport 
F e (in gyro-Bohm units) on triangularity S for elongation k = 
1.0, 1.5 and 2.0 (from top to bottom). 

In Section VI we will present results showing that the 
effect of flow shear enhancement by a weakening of the 
geodesic transfer mechanism is insignificant for DALF 
parameters. Also for these resistive electromagnetic edge 
parameters, the interchange dynamics in the normal cur- 
vature forcing is acting more as a catalyst than as a drive 
mechanism for the turbulence, with only weak balloon- 
ing of the fluctuating quantities observed. On the other 
hand, it has been established that local magnetic shear 
by flux surface deformation has a strong damping influ- 
ence on tokamak edge turbulence: comparing the case of 
a simple circular torus (with global shear only) to a case 
where the local magetic shear structure is typical of a 
shaped tokamak (but the curvature properties are those 
of a circular tor us), a reduction of transport to around a 
half was found [22] • Moreover, the effective topological 
global shear (defined, e.g., in Eq. 29 of Ref. [54j) is also 
increased by finite aspect ratio and elongation, but is not 
significantly affected by moderate triangularity (<5 < 0.6 
or so). It can thus be argued, that the major mechanism 
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for reduction of transport by elongation in the DALF case 
is magnetic shear damping, with an additional (slightly 
lesser) influence by interchange drive reduction, and an 
only insignificant zonal flow enhancement. 

We now include ion temperature gradient (1TG) driven 
dynamics by allowing for a gradient ratio rji = 2 and 
using the full six-moment version of GEM. In the ITG 
model, the flux-surface averaged turbulent transport can 
be characterised by particle transport F n = (nv x ) and by 
heat transport Qi = Qf + Q\ d with convective compo- 
nent Qf = ^Ti(fiiU x +Ti±w x ) and conductive component 
Qi d = n((^f^+fix)u x + (ni+2fij_)w x ), where u x is the 
gyro-averaged ExB velocity and w x is its FLR correction. 
Normalisation is to standard gyro-Bohm units with the 
gradient scale length set to unity, and Tj = r,/T e is the 
temperature ratio. 

When displaying heat and particle transport in relation 
to elongation and triangularity now for the ITG model, 
we find in Fig. 0a scaling of Qi ~ k -2 6 and F n ~ k~ 2 3 . 
The reduction by elongation is thus stronger than for 
the DALF model. One evident mechanism for this in- 
creased reduction in the ITG case is the weakening of 
the interchange drive of the ITG instability acting via 
the normal curvature JC y , which is slightly reduced by 
elongation in the ballooning region of the flux tube (com- 
pare Fig. [21. In contrast to DALF dyanmics, the ITG 
turbulence is strongly driven by "unfavourable" normal 
curvature, leading to a pronounced ballooning character 
of the fluctuations, as shown in Fig. [7| 

We find as an additional mechanism in the ITG case 
a significant weakening of the geodesic transfer mecha- 
nism mediated by IC X , which is discussed in the following 
section. 

Similar as in the DALF case, both heat and particle 
transport are slightly increasing with higher triangularity 
for an elliptical cross section (k = 2) in Fig. while for 
small k no influence from triangularity is observed within 
the fluctuation error bars. 

Summarising the results of this section on the influ- 
ence of flux surface shaping on turbulent transport for 
different edge parameters, we find that the main contri- 
bution to transport reduction in both the DALF and ITG 
regimes is magnetic shear damping by elongation. 

Shear flow enhancement by the reduced geodesic trans- 
fer mechanisms is found to be relatively weak in the 
DALF case, but is becoming relevant in the ITG model, 
as can bee seen from analysing the wave number and 
frequency spectra discussed in the next section. 

For fixed gradients, the local diffusivity is proportional 
to the transport and thus to the roughly inverse quadratic 
scaling with k given above. Global transport across an 
entire magnetic flux surface can in principle be obtained 
by integration of the local transport coefficients, that 
takes the variation of surface area with elongation into 
account. 
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FIG. 5: ITG: Dependence of turbulent particle transport F n 
and heat transport Qi (in gyro-Bohm units) on elongation k 
for triangularity S — 0. 
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FIG. 6: ITG: Dependence of turbulent particle transport F n 
and heat transport Qi (in gyro-Bohm units) on triangularity 
5 for elongation k = 1 and 2. 




FIG. 7: ITG: In the shaped configuration with k = 2 the 
coupling quantity h e (z) = p e {z) — 4>{z) has a reduced squared 
amplitude compared to <f>{z), indicating a weaker interchange 
ballooning character of the turbulence than in a circular torus. 
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VI. FREQUENCY AND WAVENUMBER 
SPECTRA OF FLOWS AND FIELDS 



Time traces are resolved down to l/20th of the drift 
frequency ojdw — c s /Ii, and long computational runs 
are required to account for zonal flows with lu w 0. A 
typical number of time steps is 10 5 . 

Geodesic acoustic oscillations (GAM) may be dis- 
tinctly detectable for some parameters and are ex- 
pected for circular geometry around lu/lujjw = 
(2L±/R)y/(l+n)/2. In the DALF case for n = we 
have luqam = 0.035ludw, and in the ITG case for n = 1 
we have lugam — 0.05ludw- 

In general, the GAM frequency is determined by the 
geodesic curvature term JC X and is influenced by flux sur- 
face deformation. The geodesic curvature effect is the 
action of the curvature operator JC upon the axisymmet- 
ric component {k v = 0). In local flux tube coordinates 
we find JC X ~ (1/k), and the factor 1/n in the geodesic 
curvature term does also accordingly scale the geodesic 
acoustic frequency: stronger elongation shifts the GAM 
resonance frequency in the spectrum closer to the zero- 
frequency zonal flow. The same argument, now applied 
to the normal curvature term K. v , also accounts for a n~ l 
scaling of the (interchange) drift wave frequency. 

For DALF parameters used in our computations, the 
GAM peak does not protrude very distinctly but can still 
be identified on the flat top of the spectra. In Fig. [S] the 
discrete Fourier transform spectrum (f>(u>) for a measure- 
ment of 4> at a local point in the centre of the xj/z-domain 
is shown: the spectrum below the GAM resonance is 
mostly flat except for the distinct zonal flow peak at 
lu = 0. The higher frequency part shows a typical lu~° 
cascade structure, where the exponent a changes around 
the drift frequency Cu 1, which for our choice of edge 
parameters coincides also with the Alfven frequency. 

One may identify three ranges in the spectrum: A flat 
top region between the zonal flows and GAM frequencies, 
one intermediate range between the GAM and drift wave 
frequencies with a « 1 — 2, and a cascade range between 
the drift wave and the dissipation times scale with a ~ 
3 — 5. As expected, both GAM and drift wave frequencies 
are reduced by increasing elongation to lower values lu — > 
uj/k. The whole spectrum is thus shifted to the left for 
k = 2 in Fig. 00 

Using ITG parameters, the GAM peak is more clearly 
visible for the circular torus case but vanishes for elon- 
gation k = 2, as can be seen in Fig. EI The whole spec- 
trum experiences again a shift to lower frequencies for 
the elongated case. The lu — zonal flow component in 
the ITG spectrum exhibits significantly higher amplitude 
for elongation k — 2 than for the circular torus. 
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FIG. 8: DALF: Local frequency spectra <f>(ui) for k — 1 (thin 
line) and k = 2 (bold line), both for (5 = 0. Arrows indicate 
the drift wave and GAM resonances for re = 1. 
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FIG. 9: ITG: Local frequency spectra 4>(ui) for k — 1 (thin 
line) and k = 2 (bold line), both for 8 = 0. Arrows indicate 
the drift wave and GAM resonances for re = 1. 

In addition to computing the flux surface shapes with 
varying elongation and triangularity discussed above, 
HELENA is also used to obtain equilibria with an outer 
boundary defined by an experimentally reconstructed 
separatrix position. In this way we obtain an equilib- 
rium in the typical shape of an ASDEX Upgrade (AUG) 
|23| lower single null divertor configuration (as shown in 
Fig. . The additional asymmetric shaping by the pres- 
ence of an X-point in the lower part of the torus leads to 
an increase of local magnetic shear near z — — 7r/2 and 
an associated local reduction of geodesic curvature. The 
elongation and triangularity of this configuration are oth- 
erwise comparable to the kappa-delta model with k ~ 1.6 
and S — 0.3. This AUG configuration was first used 
in edge turbulence computations in Ref. [l^, where the 
dependence of transport on j3 was established: it was 
found that the onset of MHD ballooning mode turbu- 
lence is prevented for typical tokamak edge parameters 
due to the shaping effects in realistic geometry, and the 
nature of transport is still basically of the drift- Alfven 
wave character. 
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The AUG model has, for the present flux tube position, 
properties relevant to the turbulence dynamics that are 
in some aspects in between those of up-down symmetric 
configurations with k = f.6 and 2.0. Transport in the 
AUG model is for DALF parameters found to be reduced 
to F e — 0.5 compared to F e — 0.9 for the circular torus. 
The nonlinear DALF growth rate Tm = F e /2E tot is only 
slightly lower for the AUG model with T N = 0.007±0.001 
than for the circular torus with r^r = 0.0085 ± 0.001. 

For ITG parameters, the particle transport is nearly 
identical in the AUG and SCT cases, both at F e = 1.0. 
The nonlinear ITG growth rate is also similar, r^v = 
0.0080 ± 0.0020 for AUG and Tn — 0.0088 ± 0.0024 for 
the circular torus. 

A reduction of geodesic curvature by elongation and 
X-point shaping might be expected to lead to a weak- 
ening of the geodesic transfer coupling mechanism for 
energy from zonal flows to GAMs [4!j • This can be anal- 
ysed by comparing the flux surface averaged vorticity 
(Qe(x)) — d x {v v E (x)) for the configurations of the cir- 
cular torus with the AUG model. The zonal vorticity 
(il) /rjy, set in relation to the nonlinear drive rate r^r, 
thus represents the radial shearing of zonal flows, which is 
considered responsible for the turbulent shear flow decor- 
relation and energetic damping of vortices |37| . 

Concerning the geometric effect on this, we find differ- 
ent results for the DALF and ITG parameters. 

In the DALF case, the cumulative amplitude of the 
components in the radial spectrum of flow shear (Q) /Tn 
is the same within fluctuation error bars for AUG geome- 
try as for the circular torus. The lowest k x component in 
these spectra, which is found to be stronger in the AUG 
case, makes up for only approximately 8% of the total 
vorticity. The spectra are shown in Fig. EH Flow shear 
in the DALF case thus is not significantly enhanced by 
flux surface shaping. Also the amplitude of zonal flows 
in the DALF frequency spectra in Fig. [S] is found to be 
comparable for circular and shaped configurations. The 
reduction of turbulent edge transport is in this case there- 
fore mainly a result of the local magnetic shear effect [2^] • 

In the ITG case, the flow shear rate deduced from the 
k x spectra of zonal vorticity is cleary higher in all compo- 
nents for the AUG tokamak than for the circular torus, 
as can be seen in Fig. 1111 The cumulative amplitude in 
the spectra is by 70% higher in AUG geometry than in 
simple torus geometry. Flux surface shaping introduced 
by elongation and X-point shaping clearly enhances flow 
shear significantly for ITG edge parameters. The zonal 
flow amplitude in Figgis also by a factor of three larger 
in the elongated configuration. 

This shear flow enhancement contributes to the over- 
all turbulent transport reduction by elongation, as has 
been discussed in section V, in addition to the effects of 
magnetic shear damping and the reduction of interchange 
drive. The detailed relative importance of these three 
mechanisms may however change with varying plasma 
parameters. 
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FIG. 10: DALF: k x spectra of flow shear rate in AUG geom- 
etry (bold line) and for the simple circular torus (thin line). 
The cumulative flow shear amplitude in the spectra is similar 
for both cases. 




FIG. 11: ITG: k x spectra of zonal vorticity related to non- 
linear growth rate (setting the flow shear) in AUG and SCT 
geometry. Shear flows are cleary enhanced in the AUG con- 
figuration in all components. 



VII. CONCLUSIONS AND OUTLOOK 

We have discussed the effects of flux surface shaping on 
tokamak edge turbulence and flows. Elongation strongly 
reduces the turbulent transport, whereas triangularity 
was found to have only weak (transport enhancing) in- 
fluence. Two regimes of edge plasma parameters were 
studied: an electromagnetic resistive drift wave (DALF) 
case, and the (ITG) case with an additional ion temper- 
ature gradient of r\i = 2. 

The scaling of local transport is roughly with the in- 
verse square of elongation, showing a stronger reduction 
for ITG parameters than for DALF parameters. The ITG 
result is consistent with earlier findings with the previous 
version of the GEM model [2(|. 

The mechanisms of flux surface shaping effects on tur- 
bulence have been discussed. The major cause for trans- 
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port reduction by elongation is damping of turbulence 
by stronger local magnetic shear at the upper and lower 
parts of the torus on top of increased global magnetic 
shear. The supervening reduction of the interchange 
drive by a lower normal curvature has an additional tur- 
bulence reduction effect in particular for ITG parameters, 
where the interchange drive is slightly weakend. In the 
ITG case a significant enhancement of zonal flow shear 
by elongation and X-point shaping further contributes to 
transport reduction. The shear flow effect is on the other 
hand found to be negligible for DALF parameters. 

Both the DALF and ITG parameter regimes used in 
our computations seem consistent with the presently still 
poor experimental knowledge on edge ion temperature 
profiles, though this situation is rapidly improving |55| . 
It can be assumed that additional core plasma heating, 
e.g. by neutral beam injection (NBI) used in particu- 
lar for obtaining the transition to an H-mode, is both 
increasing the ion temperature at the pedestal top and 
the ion temperature gradient in the overall edge pedestal 
region. Based on our results, the shaping of tokamak 
flux surfaces by elongation and by the additional pres- 
ence of an X-point like in AUG is expected to facilitate 
an enhancement of zonal flows when the role of the ion 
temperature gradient dynamics is pronounced by addi- 
tional heating. 

It can be expected that this enhancement of zonal flows 
found for ITG parameters is even more pronounced for 



more strongly shaped flux surfaces nearer to the separa- 
trix. Preliminary results suggest that such initial radially 
local flows have the ability to spread over wider regions of 
the pedestal in simulations with global profile evolution. 

Some L-H transition theories include zonal flows as a 
trigger mechanism to induce the transition to a sustained 
mean shear flow. The enhancement of zonal flows by 
flux surface and X-point shaping in the presence of ion 
temperature gradient steepening may offer an additional 
explanation for the observation that H-mode access in 
tokamak experiments is facilitated by the presence of a 
divertor. 

The possibility to obtain a confinement transition 
within first principle computations of edge turbulence 
will thus have to be studied with a code that includes 
full temperature dynamics, realistic flux surface geome- 
try, global profile evolution, and a coupling of edge and 
SOL regions including realistic sheath boundary condi- 
tions, while in addition it maintains sufficient grid reso- 
lution, grid deformation mitigation, and energy plus en- 
strophy conservation in the vortex/flow system. 
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